


#Exercicio 9.2 - Analise exploratoria simultanea de duas variaveis

install.packages("plotrix")
library("plotrix")

#dados para teste:
x=rnorm(30,8,2)
y=rnorm(30,7,2)

a.expl=function(x,y)
{
  if 
  (length(x)!=length((y)))
    {
    stop("Vetores de tamanhos diferentes")
    }
  x11()
  par(mfrow=c(2,2))
  library("plotrix")
  z = list(x,y)
  multhist(z, xlab= "Variavel x", ylab="Variavel y", col=c(1,8) )     #histograma com as 2 variaveis
  boxplot(x,y, xlab= "Variavel x", ylab="Variavel y")    
  axis(side=1, at=1:2, labels=c("x", "y"))   #legenda
  qqnorm(y, col="red",main="Grafico Quantil x Quantil",xlab="Quantil Teorico",ylab="Quantil da Amostra")
  qqline(y)
  mod = lm (y~x)   #modelo linear
  plot(y~x)
  abline(mod, col=2)
  par(mfrow=c(1,1))
  summ.x=summary(x)
  summ.y=summary(y)
  coeff.int= coef(mod) [[1]]    #isolando o intercepto (1)
  coeff.inc= coef(mod)[[2]]     #isolando a inclinacao da reta (2)
  resultado=list(summ.x, summ.y, coeff.int, coeff.inc)
  names(resultado)=(paste(list("Sumario da variavel x", "Sumario da variavel y", "Intercepto", "Inclinacao da reta")))
  return(resultado)
}

a.expl(x,y)
